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We investigated the noise of interferometric gravitational wave detectors due to heat energy de- 
posited by cosmic-ray particles. We derived a general formula that describes the response of a mirror 
against a cosmic-ray passage. We found that there are differences in the comic-ray responses (the 
dependence of temperature and cosmic-ray track position) in cases of interferometric and resonant 
gravitational wave detectors. The power spectral density of vibrations caused by low-energy sec- 
ondary muons is 100-times smaller than the goal sensitivity of future second-generation interferom- 
eter projects, such as LCGT and Advanced LIGO. The arrival frequency of high-energy cosmic-ray 
muons that generate enough large showers inside mirrors of LCGT and Advanced LIGO is one per 
a millennium. We also discuss the probability of exotic-particle detection with interferometers. 

PACS numbers: 04.80.Nn, 95.35.+d, 95.55. Vj, 95.85.Ry 



I. INTRODUCTION 



Recent improvements of the sensitivity and operational 
stability of gravitational wave detectors is remarkable. 
Observation runs have already been performed in sev- 
eral interferometer (LIGO [1], VIRGO 0, GEO 0, 
TAMA 0, CLIO 0) and resonator (ALLEGRO Q, EX- 
PLORER 0, NAUTILUS0, AURIGA 0, NIOBE \v^, 
MARIO SCHENBERG [nf) projects. In order of gravi- 
tational wave detection, the reduction of noise and fake 
triggers is crucial, since the gravitational wave ampli- 
tude and number of events are expected to be small and 
rare. In 1969, it was pointed out that cosmic-ray particles 
could cause fake triggers in resonators 12]. The interpre- 
tation of excitations of resonators by cosmic-ray particles 
is follows: The heat energy deposited by cosmic-ray par- 
ticle passages induces temperature gradients around their 
tracks, and thermal stress excites internal vibrations of 
the resonator. These phenomena have been investigated, 
for example, observations of excited resonator vibrations 
by beams from accelerators [HI, HH, simultaneous de- 
tection of resonator excitation and cosmic-ray particles 
[TH , and studies of exotic events in super-conductive res- 
onators [H, [I?], [H| ■ In some research, resonators were 
also operated and treated as exotic-particle detectors 
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These studies suggest that cosmic-ray 
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heating is a possible noise source in interferometric de- 
tectors [23|, 0, US HH (other effects on interferometers, 
the momentum and electrical charge brought by cosmic- 
ray particles, are discussed in Refs. [23l 124 125I |27| ) . 

We investigated details of this effect by cosmic-ray en- 
ergy deposition in interferometers. A formula that de- 
scribes the response of a mirror against a cosmic-ray 
passage was derived. This formula reveals differences 
between the cosmic-ray responses of interferometers and 
resonators. We used it to evaluate the amplitude of vi- 
brations caused by cosmic-ray particles in typical cases of 
interferometers, and examined the effect in gravitational 
wave detection. We also considered the probability of 
exotic-particle detection with interferometers. 



II. FORMULA OF EXCITED MOTION BY A 
COSMIC-RAY PARTICLE 

A. Outline of derivation of the formula 

In order to simplify the discussion, the mirror vibra- 
tion excited by a cosmic-ray particle is investigated. A 
vibration excited by many particles, like a shower, is a 
superposition of that by one particle. The excitation by 
a particle is considered under the following assumptions. 
The particle goes straight and never stops in the mirror. 
Its speed is faster than that of sound in the mirror. A 
long and narrow heated volume appears at the instant of 
particle passage. 

The heat-conduction equation is solved in order to cal- 
culate the time evolution of the temperature gradient. 
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The vibration of the mirror is examined using the equa- 
tion of motion of an elastic body with thermal stress, 
which is proportional to the thermal gradient. 



B. Formula 

Since the heated volume is smaller than that of the 
mirror, itself [28l . [29| . the mirror and the initial heated 
volume are treated as an infinite body and a line, respec- 
tively. The direction of the cosmic-ray track is taken as 
the z-axis. The heat-conduction equation is described as 

M 



d_ 

di 



ST 



P C 



A(6T) 



1 dE 
~p~C~dJ 



S(x)S(y)S(t), (1) 



where ST is the temperature difference caused by a 
cosmic-ray particle. The quantities K,p,C and dE/dl 
are the thermal conductivity, density, specific heat per 
unit mass, and energy loss of a particle per unit length, 
respectively. The solution is described as [3l| 
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The radius of the heated volume increases with time due 
to conduction. The time when the heated area radius 
becomes a is 



pCa 2 
Ak 



(3) 



The equation of motion of an elastic body with thermal 
stress is described as I30B 



d 2 u 
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2(l + cr)(l-2cr) 
Ya 



1 -2a 



grad divit 
div<5T, (4) 



where u represents the displacement of a volume ele- 
ment in the elastic body. The quantities Y, a and a 
are Young's modulus, the Poisson ratio, and the linear 
thermal-expansion coefficient, respectively. By substi- 
tuting Eq. ([2]) for Eq. (f4]), we obtain the output of a 
interferometer X, 



V = / u opt {r)P(r)dS, 

i surface 



(5) 



where it op t is the optical axis component of u and P is 
the intensity profile of the laser beam, 



P(r) 



7T7V 



■ exp 



2r2 
r Q 2 



(6) 



The quantities r and r$ are the distance from the optical 
axis and the beam radius. We employ the modal expan- 
sion method [HI, HH, [34| to calculate u and X. In this 



method, u and X are represented by a superposition of 
the resonant mode displacement, 



n 



(7) 
(8) 



where w n and q n represent the displacement and time 
development of the n-th resonant mode, respectively. 
These basis functions are normalized to satisfy a con- 
dition [HE EH, 



/ w n , opt (r)P(r)dS = 1, 

J surface 



(9) 



where u>„ i0 pt is the optical axis component of w n . The 
equation of motion of each mode is the same as that of a 
harmonic oscillator, 

-m n u 2 q n (uj)+m n uj n 2 [l + i(fi n (uj)}q n (uj) = F n (uj), (10) 

in the frequency domain. The quantity <p n is the loss an- 
gle, which represents dissipation of the n-th mode [Hj]. 
The force F n applied on the n-th mode is related to the 
thermal stress. The quantities m n and w n are the effec- 
tive mass and the resonant angular frequency [33], [33 |. 
The effective mass is defined as 



pw n (r) ■ w n (r)dV. 



(11) 



volume 



The quantities q n {u) and F n (u>) are the Fourier compo- 
nents of q n and F n , respectively: 



X{t) 



1 

2^ 



X(t)exp(-iojt)dt, (12) 
X(to) exp(icot)duj . (13) 



The force F n is obtained from the modal decomposi- 
tion of the thermal stress on the right-hand side of Eq. 
(0]). The decomposition procedure [33|, H3 is as follows. 
The thermal-stress term is multiplied by w n . The inte- 
gral of this inner product over all the volume is F n . This 
force F n decreases after the heated volume scale, a, be- 
comes larger than the n-th mode wavelength. In order 
to simplify the discussion, it is assumed that the time 
evolution of F n [l^, [35[ is expressed as 

(t) = U(0) exp (-£) (t>0) 
{ (t < 0) 

The quantities F n (0) and r„ are the initial value and the 
decay time of the force, respectively. The initial value, 
F n (0), is written as [13, HI 



F n (0) 



Ya 1 
l-2cr^C 



. dE , , 
divw n dl ) —. (15) 
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The integral along the cosmic-ray track represents the 
coupling between the thermal stress and the n-th mode. 
The coefficient 1 / (pC) is a factor used to transform the 
heat energy into the temperature gradient. The force F n 
is described as a product of the temperature gradient and 
Ya/(1 — 2er). From Eq. ([3|), the time r„ when the heated 
volume radius becomes comparable to the wavelength of 
the n-th mode is expressed as 



pCX n 
An 



ir 2 P Cv 2 ir 2 YC 



(16) 



The quantities A„ and v are the wavelength and sound 
velocity: 
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The Fourier component of F n in Eq. (] 14[) is written in 
the form 



Fn{0) Tn 

2ir 1 + \ujT n 



(19) 



We now write down the formula of the mirror vibration 
excited by a cosmic-ray particle using Eqs. ([8]), (flUf . (fT5|) 
and IT911: 



F n (u) 



n 

= E 



m n uj 2 + m n uj n 2 (l + i<j) n ) 
1 Ya 1 dE 



2vr 1 - 2cr pC dl 

*E 



i 



-m„uj 2 + m n u n 2 (l + i(j) n ) 1 + \LOT n 
\ j / divw n dl ) . (20) 



C. Frequency dependence of the formula 

A schematic view of the frequency dependence of the 
modes q n (oj) in Eq. (I20|) is shown in Fig. Q] Here, 
we discuss the frequency dependence below the resonant 
frequencies of the mirrors, because the observation band 
of interferometers (around 100 Hz) is below the funda- 
mental mode (the order of 10 kHz). The cut-off fre- 
quency, 1/(27tt„), is extremely smaller than the resonant 
frequency, uj n /{2ir), as shown in Fig. [TJ because sound 
is generally faster than heat conduction. The absolute 
value |g„(o;)| is inversely proportional to the frequency 
between 1/(27tt„) and w„/(27r). Below the cut-off fre- 
quency, 1/(27tt„), |<7„(w)| is constant. 



FIG. 1: A schematic view of the frequency dependence of 
the Fourier components of mode motion excited by a cosmic- 
ray particle, q n (u>) in Eq. I|20|) . The absolute value |q„(u;)| is 
inversely proportional to the frequency between l/(27rr n ) and 
u>„/(2-k). Below the cut-off frequency, l/(27rr n ), it is constant. 
The "highest" mode is that with a wavelength comparable 
to the beam radius, ro, and contributions of higher modes 
are negligible in the summation of Eq. (|20[) 37, 38]. The 



cut-off frequency of the "highest" mode, 1/(27tt,- ), is smaller 
than the fundamental mode resonant frequency, uii/(2n), in 
general. The cut-off frequency, 1/(2ttt„), of the lower mode 
is smaller. In the range between l/(27r7v ) and u>i/(2tt) (the 
"high" frequency region in this graph), |X(w)| = |X]<?n(w)| 
is inversely proportional to the frequency. Below the cut- 
off frequency of the fundamental mode, 1/(2ttti) (the "low" 
frequency region in this graph), |X(oj)| is independent of the 
frequency. The approximation formulae in the "high" and 
"low" frequency regions, Eqs. (|30|l and (|34p . are derived using 
Eq. J20J. 



The frequency dependence of = J2q n (u;) is as 

follows. The "highest" mode in Fig. Q] is that with a 
wavelength comparable to the beam radius, ro, and con- 
tributions of higher modes are negligible in the summa- 
tion of Eq. (J20J) [13, HI]. The thermal relaxation time, 
r ro , for this " highest" mode is described as 



pCrl 
An 



(21) 



from Eq. (fl6|) . The cut-off frequency of the "highest" 
mode, l/(27rr ro ), is smaller than the fundamental mode 
resonant frequency, aq/(27r), in general, as shown in Fig. 
[U The cut-off frequency, 1/(27tt„), of the lower mode is 
smaller from Eq. (fTo) . In the range between l/(27rr ro ) 
and uji/(2tt) (the "high" frequency region in Fig. [1]), 

l^( w )l — IS"?"^)! i s inversely proportional to the fre- 
quency. Below the cut-off frequency of the fundamental 
mode, 1/(27tti) (the "low" frequency region in Fig. [1]), 
|X(a>)| is independent of the frequency. From Eq. (JTC 
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the relaxation time, n, is described as 



Tl 



pCXx 2 _ pCR 2 
An K 



(22) 



because the wavelength of the fundamental mode, Ai, is 
comparable to the mirror diameter, 2R. 

The typical values of \X(u>)\ in the "high" and "low" 
frequency regions of Fig. Q] are evaluated using Eq. 
([2"0]). In the "high" frequency region, l/(27rr ro ) < / < 
u>i/(2tt), it is approximated as (\(f> n \ -C 1 in usual cases) 



\X(u)\ 



E 



1 Ya 1 dE 1 
2n 1 — 2(7 pC a-> ^ — ' m n w n ' 

l/(27TT ro ) < /< Wl /(27T). 



divu>„ e?Z 



(23) 



The sign of the integral in Eq. (|23[) depends on the 
modes. The typical absolute value of the summation in 
Eq. (f23|) is evaluated as the square root of a summation 
of squares of the terms. Equation (|23|) is rewritten as 



1 Fa 1 rfE 1 1 
2tt 1 - 2a~pC~dlui 



1 



m„ 2 w„ 4 



diviiL, cM 



(24) 



l/(27TT ro ) < / < Wi/(27T). 



The integral along the cosmic-ray track in Eq. (|24|) is 
evaluated as follows. The average of the length of the 
cosmic-ray track is comparable to the radius of the mir- 
ror, R. The average of |io n | 2 , (l^nl 2 ), is related to Eq. 

CD. 

m„ = J P \w n \ 2 dV = M (\w n \ 2 ) , (25) 
where M is the mass of the mirror. This equation gives 

V<RP) = ^. (26) 

The divergence of w n in Eq. (f24|) can be represented by 
the product of w n and the wavenumber Lo n /v, because 
w n is the basis of the solution of the wave equation. Con- 
sequently, Eq. ([24l is described as 



\X(u)\ 



1 Ya 1 dE R 1 
2tt 1 - 2a pC dl V \J~M w ' 

l/(27TT ro ) < / < Wl/(27T). 



(27) 



The summation in Eq. (|27p is the same as the response 
of a mirror against a static force [39|, |40| , 
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m n Lu n 2 y/nYro 
it 

In order to simplify the discussion, the relation 

M = irpR 3 (29) 



is assumed. The radius of the mirror, R, is nearly equal 
to its thickness, H, in usual cases of interferometric grav- 
itational wave detectors. Using Eqs. (fl8|). (|28|) and ([29jl . 
Eq. (|2T|) is written in the form 



1 ay/1 - a 2 I dE 1 I 
2tt 7 /4 l-2cr ^~dl JRF^ 
l/(27rr ro ) < / < wi/(2tt). (30) 



Equation (|20|) in the "low" frequency band of Fig. [T] 
is evaluated in the same manner as in the previous para- 
graph. Using Eq. (|T6|) . the result is written as 



\X(u>)\ 



1 Ya 1 dE 
2tt 1 - 2a~p~C~dT 



\ n v 

V7r^ 3/2 a 1 d£ 1 

~2~ 1 - 2a^~dT~^/RV 2-j n 

f < l/(27rn). 



(31) 



It can be seen that only the fundamental mode is domi- 
nant, because of the frequency dependence of w„ 6 . The 
quantities of this mode are as follows [2(| Hi" 



U3\ 




(32) 
(33) 



Equation ([31]) is expressed as 
1 



\X{u>)\ 



v^tt 3 1 - 2a k dl 



1 rl 

--R, f<l/(2nn). (34) 



D. Formula in time domain 

Here, we discuss the excitation formula Eq. (|2"0"|) in 
the time domain. Only a contribution of the n-th mode 
is considered in order to simplify the discussion. If the 
Q- value, Q n — l/<p n (u> n ), is larger than unity and l/r„ is 
smaller than u) n , q n in the time domain (t > 0) is written 
in a form 



,,x F n (0) ( t 

q n {t) ~ j ex P — 



F»(0) 
m n u) n 2 



cos(w„i) exp 



2Qn 



(35) 



Figure shows the fundamental mode, qi, in the time 
domain. 

The second term in Eq. (|35|) is dominated by X(ui) 
near the resonant frequency. This is the excited resonant 
vibration and its decay. The outputs of resonant detec- 
tors are described with this term. The first term in Eq. 
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FIG. 2: Vibration of the fundamental mode caused by a low- 
energy cosmic-ray muon in the time domain, qi in Eq. (I35|l . 
In the calculation, the material values of sapphire at room 
temperature are used. Since the Q-value is extremely high, 
i.e. the decay time is longer than the period of the resonant 
motion, we are not able to see the resonant motion of the one 
period in this graph. 



(|35|) represents the drift of the center of the resonant vi- 
bration caused by the relaxation of thermal stress. This is 
dominated by X(oj) below the fundamental mode. The 
outputs of interferometric detectors are described with 
this term. 

The initial amplitude of ^ n 2F n (0)/(m„w„ 2 ) is evalu- 
ated as Eq. ([30]) . 
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1 dE 
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(36) 



Here, we consider a sapphire mirror at room temperature. 



III. DISCUSSION ABOUT THE FORMULA 
A. Effect on low-temperature interferometers 

In some future projects using interferometric detectors 
as LCGT I4l| and ET [431 , mirrors will be cooled in order 



to reduce the thermal noise (for example, LCGT mirrors 
at 20 K). In the quantities of the force, F n , which is 
related to the thermal stress, a, C and k in Eqs. (fT5|) 
and (Unj) strongly depend on the temperature [43]. The 
initial value of the force, F„(0) in Eq. (fT5|) , is propor- 
tional to a/C. The decay time of the force, r„ in Eq . 
((THJ), is proportional to C/k. The Griineisen relation [44[ 
predicts that the ratio a/C is independent of the temper- 
ature. The initial force, F n (0), and the initial amplitude 
of the excited vibration do not depend on temperature. 
On the contrary, in the case of crystals, the decay time, 
t„ (oc C/k), of the first term in Eq. (|35[) is extremely 
short at the cryogenic temperature, because of the small 
C and large k. The cut-off frequency, l/(27rr„), in the 
low-temperature region is higher than that at room tem- 
perature (e.g. Ref. [111). For example, the cut-off fre- 
quencies of sapphire at room temperature, obtained from 
Eqs. flUD and (JH]), are: 



27TT1 



0.13 mHz 



25 cm N 



4 g/cm 3 \ /7.9 x 10 2 J/kg/K 



C 



2R 



40 W/m/K 



(37) 



27T7V 



= 9.0 mHz 



3 cm* 



4 g/cm 3 \ ^7.9 x 10 2 J/kg/K 
C 



r J \ 40 W/m/K 
The values at 20 K are: 



(38) 



1 



27TT1 



= 58 Hz 



4 g/cm 3 \ /0.69 J/kg/K 



25 cm 



C 



2irT r , 



= 4.0 kHz 



3 cm 



2R J \ 1.6 x 10 4 W/m/K 
4g/cm 3 \ /0.69 J/kg/K 



(39) 



C 



J'o 



1.6 x 10 4 W/m/K 



(40) 



At low temperature, the cut-off frequencies are near the 
observation band of gravitational wave detectors (around 
100 Hz). The "high" frequency approximation of |-X"(o>)| 
in Eq. (|30[1 is only appropriate for room-temperature 
interferometers, and not valid for cryogenic interferom- 
eters. In order to show the effect of the cooling mir- 
rors, \q\ (uj) | of a sapphire mirror at 300 K and 20 K are 
plotted in Fig. [3] In the low-frequency region, X(uj) be- 
comes much smaller due to cooling. The "low" frequency 
approximation of in Eq. Q34p is proportional to 

ol/k = a/C x C/k oc C/k, which is small in the low- 
temperature region. Since the decay time of the force F n 
becomes shorter, it is difficult to excite the low- frequency 
component. The Fourier components |gi(w)| in the high- 
frequency region of Fig. [3] are comparable in the cases 
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FIG. 3: Fourier components of the fundamental mode motion, 
|gi(w)t, of a sapphire mirror excited by a low-energy cosmic- 
ray muon at 300 K (solid line) and 20 K (dashed line). The 
cut-off frequency, 1/(2ttt\), at 20 K is higher than that at 
300 K. The mirror cooling reduces the low-frequency compo- 
nent because the decay time of F\, which is related to the 
thermal stress, is shorter. The higher frequency component is 
independent of temperature because of the Griineisen relation 

H. 




A>B 



Bar resonator 




1 A=B 



FIG. 4: Track position dependence of the cosmic-ray heating 
effect of interferometers and bar resonators. In the outputs of 
interferometers, the vibration caused by a particle that goes 
along track A near the beam spot is larger than that along 
track B far from the beam spot. On the contrary, in the case 
of bars, the vibration by track A is the same as that by track B 
because the displacement of the resonant modes is symmetric 
or antisymmetric with respect to the center of the bar. The 
dashed line shows the fundamental mode deformation. 



of 20 K and 300 K. The "high" frequency approximation 
of |X(w)| in Eq. (|30|) is independent of temperature be- 
cause it is proportional to a/C and related to the first 
term of Eq. (f35|) at t ~ 0. From the cut-off frequencies in 
Eqs. (f3T)]) and (f^D|) and a comparison between Eq. (IM)) 
at 20 K and Eq. (|30|) . it can be seen that the vibra- 
tion of a cooled sapphire mirror excited by a cosmic-ray 
particle in the observation band (around 100 Hz) is a 
few-times smaller than that at room temperature. This 
is an advantage of cryogenic interferometers in addition 
to the suppression of the thermal noise [HI, [46|, IttI. t her- 
mal lensing effect [48[ and parametric instability [4{J. It 
must be noted that motions excited by cosmic-ray par- 
ticles in resonant detectors are independent of tempera- 
ture [l6|, [I?], EH- This is because the initial amplitude, 
the second term of Eq. (|3"5|) at t ~ 0, does not depend 
on temperature. 



B. Cosmic-ray track position dependence 

In the calculation of Eq. (|2"D|) , the signs of the inte- 
gral terms in the summation are important. The sign 
depends on the positions of a cosmic-ray track and the 
laser beam spot, because the displacement of the mode 
w n is normalized to satisfy Eq. © |50(. If the particle 
track is near the beam spot, the signs of the integrals of 
many modes are the same, because the basis functions, 
w ni on the track are similar. If the track is far from the 
spot, w n on the track and the integral signs are different 



for various modes. Since the sign of q n in Eq. ([20]) below 
the fundamental mode is the same as that of the inte- 
gral, |A(w)| below the first mode is larger and smaller if 
the cosmic-ray track is near and far from the beam spot, 
respectively. Another explanation about the cosmic-ray 
track position dependence is as follows. The heated vol- 
ume on the particle track pushes around them. Since 
the center of the mirror does not move because of the 
conservation of momentum, a larger motion is observed 
if the track is near the beam spot. The track position 
dependence of the cosmic-ray heating effect in interfer- 
ometers is different from that in bar resonators, as shown 
in Fig. 2J In the outputs of bars, the vibration caused 
by a particle that goes along track A in Fig. 2] is the 
same as that along track B, because the displacement of 
the resonant modes is symmetric or antisymmetric with 
respect to the center of the bar (the dashed line in Fig. 
U shows the fundamental mode deformation) . The dif- 
ference between interferometers and the bar resonators 
is related to the number of modes to be considered. In 
the case of bars, only the fundamental mode is taken into 
account. On the other hand, in the case of interferome- 
ters, many modes contribute to the response of a mirror. 
The signs of these modes have an important role. The 
discussion above is the same as that about thermal noise 
below the fundamental mode caused by inhomogeneously 
distributed loss [H, [H, HH . 
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IV. APPLICATION I — LOW-ENERGY 
COSMIC-RAY PARTICLES 



Low-energy cosmic-ray particles and interaction 
with matter 



Primary cosmic rays generate extensive air showers 
in the atmosphere. Cosmic-ray particles on the ground 
are secondaries from air showers. Three quarters of sec- 
ondary particles at sea-level are muons. The remainder 
are almost electrons (52j . Muons with small energy (less 
than about 0.22 GeV) and electrons can be neglected be- 
cause it is difficult to penetrate matter around the mir- 
rors, for example, walls of buildings, vacuum chambers 
[521 ]. The speed of muons that arrive at the mirrors is 
comparable to that of light. The flux of these cosmic-ray 
muons at sea-level is about 2 x 1CP 2 /cm 2 /sec [H, 

Since the number of higher energy muons is smaller 
[54j , the energy of most of the cosmic-ray muons is be- 
low 100 GeV. In this low-energy region, the dissipation 
process in material is dominated by ionization (54 . |55| , 
which is Coulomb scattering with electrons in atoms of 
matter [591 . The ionization loss is about 



= 2 MeV/(g cm" 2 ), 

p at 



(41) 



and almost independent of the particle energies [5J, [55j . 
The typical loss per unit length, dE/dl, is several 
MeV/cm. 

The effect of mirror excitation by cosmic-ray particles 
depends on the arrival frequency of particles and the de- 
cay time of the vibrations. If the decay time is longer 
than the interval of the particle arrivals, the mirror vi- 
bration is maintained. If the next muon comes after the 
vibration has disappeared, the vibration can be treated 
as a burst event. The number of muons, N, that hit a 
mirror at sea-level per unit time is expressed as 



N = 2 x 10 -2 /cm 2 /sec x 2R x H 
2R \ ( H 



= 8 /sec 



25 cm 7 V 15 cm 



The average arrival interval of muons, l/N, is 



1 /25 cm\ / 15 cm 

N = °- 13sec [-wr )hr 



(42) 



(43) 



The decay time of the fundamental resonant vibration is 
described as 



7r/l 



x 10 sec 



40 kHz 



) (9l 

h A10 8 



(44) 



Since the Q-values of mirrors used for gravitational wave 
detectors are at least 10 6 , the decay time is extremely 
larger than the expected arrival interval of cosmic-ray 
particles. 



B. Power spectral density 

The power spectral density, G cos (/), of vibrations 
caused by low-energy cosmic-ray particles has been cal- 
culated (e.g. fiefs. [H, Hi|). It is assumed that arrival 
time of particles and track position in a mirror are at 
random. Since there are four mirrors in an interferome- 
ter, the one-side power spectral density of the noise of an 
interferometer output is written in the form [56] 



32it 2 N 
L 2 



(45) 

where L is the length of the interferometer arms. The 

quantity ^|X(oj)| 2 ^) is the ensemble average of |X(o>)| 2 , 

which is the vibration caused by a muon. To evaluate 
the power spectrum of room-temperature interferome- 
ters, the square of Eq. (|3TJ)) is used as the ensemble aver- 
age, because this formula is appropriate to calculate the 
typical |X(w)| at 300 K and around 100 Hz, as shown 
in Sec. IhTaI From Eqs. (gTJ), © and 69, the power 
spectrum of room-temperature sapphire at sea-level is 
written as 153 



3 km 



y/Gcosif) = 1.3 x 10 




L J \5 x lO- 6 /K 
7.9 x 10 2 J/kg/K 



1/2 



C 

( 100 Hz 



\25 cm ) \ r o J V / 

The sensitivity of future second-generation interferome- 
ter projects, such as LCGT [4l| and Advanced LIGO [HI, 
is on the order of 10~ 24 /VHi at 100 Hz. Therefore, the 
effect of low-energy cosmic-ray particles is not a serious 
problem, even in these future projects. 



(46) 



APPLICATION II 



SHOWER 



High-energy cosmic-ray particles often generate many 
particles (showers). From Eqs. flJT} and (|6"3")) . if 1000 
shower particles pass in a mirror at the same time, the 
excited vibration is large enough to be detected by fu- 
ture second-generation interferometers, such as LCGT 
[4l| and Advanced LIGO [H| (e.g. Refs. 0, H H). 
Such excitations caused by cosmic-ray showers have been 
observed in a resonator [15j ■ 

We investigated the effect of a shower generated by 
a high-energy muon inside a mirror with a Monte-Carlo 
technique [59j . It was assumed that the material is sap- 
phire. We evaluated the probability that a high-energy 
muon that runs in a 30 cm thickness sapphire generates 
more than 1000 electrons. In this simulation, the flux of 
muons at sea-level was expressed as [6l| 

,-6 ,^,„J E V 2 - 7 

1 TeV/ 



7 M (> E) = 1.1 x 10~ 6 /cm 2 /sec 



(47) 
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if E is more than 1 TeV. Our simulation showed that the 
number per unit time and per a mirror of muons that 
generate more than 1000 electrons N(> lOOOe) is (the 
typical energy of such muons is about 10 TeV [6a ]) 



N(> lOOOe) = 1.0 x 10~ n /sec 



2R 
25 cm 



H 



15 cm, 

(48) 

It must be noted that this average arrival number, 
N(> 1000c), was overestimated, because only a part of 
muons has more than a 30 cm length track in a sapphire 
mirror. Since there are four mirrors in an interferometer, 
the average arrival interval is 



1 



4iV(> lOOOe) 



= 7.8 x 10 2 year 



25 cm 
2R 



15 cm 
H 



(49) 

The effect of showers generated by high-energy muons 
inside mirrors is not a serious problem. 

In the case of a shower that occurs near a mirror, the 
energy of an original particle that generates 1000 parti- 
cles is about 1 TeV [25|. Since the spread of particles 
in a TeV energy shower is quite large, the typical size 
mirror of interferometers can not contain all of the en- 
ergy of a thousand particles. In order to know how often 
more than 1000 particles go into a mirror, accurate sim- 
ulations about shower generation in apparatus around 
mirrors (for example, vacuum chambers and vibration 
isolation systems) and the response of a mirror are nec- 
essary as resonators [i^, HI] . This is our future work. 



A. Signal-to-Noise ratio of interferometers 

Since the time evolution of an excited motion by an 
exotic particle is predicted from Eq. ([20]). the matched 
filtering method can be applied to the outputs of detec- 
tors. The output of a matched filter, the signal-to-noise 
ratio (S/N), is defined as [7? 



S/N 




ISHI 2 

GdctU) 



df, 



(50) 



where S(u>) and Gdet(/) are the Fourier components of 
the signal and the one-side power spectral density of the 
noise of gravitational wave detectors, respectively. 

In the case of interferometers, G<jet(/) and S(ui) of Eq. 
([50)1 are the strain noise, Gj nt (/), and the ratio of the 
Fourier component of the motion excited by an exotic 
particle, X{ui), to the arm length, L, respectively. It is 
supposed that the temperature is 300 K. Here, we recall 
Eq. (JSPJ in the form 



A 



A vl 



1 



7r 3 / 4 Yy/Rl^ /' 

1 Ya 1 dE 



4tt 2 1 - 2a pC dl' 
S/N is expressed by using Eqs. ([50]) and ([STjl : 



S/N int = 



4 7 r 1 /4 y 4yr — ^2 



df 



/ 2 Gint(/)' 



(51) 
(52) 

(53) 



VI. APPLICATION III — EXOTIC-PARTICLE 
SEARCH 

The effect of cosmic-ray particles on gravitational wave 
detectors suggests that the detectors are useful to search 
for exotic particles that dissipate a large amount of en- 
ergy in material. Ideas that resonators can be used as 
magnetic monopole fl9| or mirror dust particle [22| de- 
tectors were proposed. The upper limits of the flux of 
nuclearite l69L l7fJl from the operation of resonators were 
reported [20|, |2l|. Here, we discuss interferometers as 
exotic-particle detectors in comparison with resonators 
(bars [1,0, MM E3)- 

In order to detect exotic particles or other rare events, 
a larger aperture and higher sensitivity are required for 
detectors. The cross section of a bar resonator is 10-times 
larger than that of an interferometer [7l| . The area of a 
bar is about 1.8 m 2 (diameter, 0.6 m; length, 3 m) [2l| . 
The cross section of four mirrors of an interferometer is 
about 0.15 m 2 (diameter, 0.25 m; thickness, 0.15 m). We 
discuss the sensitivity of interferometers and bars for an 
exotic particle passage. 



B. S/N of bar resonators 

In the case of bar detectors, S(cu) and Gdct(/) of Eq. 
([50]) are the force applied by an exotic particle to the n-th 
mode, F n (oj), in Eq. ([T5]) . and the tidal force, which cor- 
responds to the strain noise, Gbar(/), respectively. Equa- 
tion (fl"9]) is rewritten as 

\S(oj)\ = |F„(w)| - Aj divw n dl\ , (54) 

because the cut-off frequency, l/(27rr„), is lower than the 
resonant frequency. Here, we take the fundamental mode 
of bars into account. Under the same approximation in 
the derivation of Eqs. ([3"0"]) and ([M]) , S(oj) is expressed 
as 

|S( W )| = |AM|~4-~. (55) 

The tidal force that corresponds to the strain noise, 
G b ar(/), is obtained from Refs. [Illzllzli 

G det (/) = (M h u 1(h) 2 ^) G bar (/), (56) 
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where Mb, wi(b) an d I are the mass, angular resonant 
frequency of the fundamental mode and length of a bar. 
S/N is given by Eqs. ([SO ]) . ([55 ]) . and (j56]) : 



S/N bar = 




MjWi( b ) 



(57) 



C. Comparison between interferometers and bar 
resonators 

Here, we discuss the effects of an exotic particle on 
interferometers and bar detectors by using Eqs. (|53j) and 
(|57[) . The integral term only depends on the sensitivity of 
gravitational wave detectors. It must be noted that the 
weight, l// 2 , originates from the frequency dependence 
of S(lu), i.e. Eqs. (ST) and (SI). The integral term in 
Eq. |53|) for future second-generation interferometers, 
e.g. the LCGT project [U, is 



df 



/ 2 G in t(/) 



= 3.0 x 10 2 



(58) 



The typical goal sensitivity of bar resonators [73] is 
\J Gbar(/) ~3x 10~ 22 /vTfz in the frequency range be- 
tween 850 Hz and 950 Hz. The integral term in Eq. ([57) 



4f 



/ 2 Gbar(/) 



3.7 x 10 



1!) 



(59) 



The integral term of interferometers is 1000-times larger, 
because interferometers have higher sensitivity and a 
wider observation band. Since the observation band of in- 
terferometers is lower than that of resonators, the weight- 
ing function, l// 2 , increases the integral term of inter- 
ferometers. 

The factors, except for the integral term and A in Eqs. 
(|53|) and (|57|) , are the ratios of the responses to an exotic 
particle to that to the gravitational wave. If this factor 
is large, the detector is more suitable for exotic-particle 
searches. This factor of interferometers is 17811 



- = 4.9 x 10~ 14 /N ' ' 



0.96 / 
4 x 10 11 Pa\ /25 cm\ 1/2 



Y 



2R 



6 cm\ 1 ^ 2 / 3 km 



In the case of bar resonators, this factor is 
2x/2tt 4 



(60) 



M)Wl(b) 2 ^ 



1.2 x 10 J /N 



2300 kg 

Mb 



900 x 2tt rad/Hz\ /3 m 



^i(b) 



— ).(6D 



The factor of bars is extremely larger. The main reason 
for this difference comes from the sizes of the detectors, 
L and I, An exotic-particle detector must be a good 
displacement sensor. A smaller size detector is a bet- 
ter displacement sensor, if the strain (gravitational wave) 
sensitivity is the same. The better strain sensitivity of 
interferometers shown by Eqs. ([58]) and (l59|) is canceled 
by their larger size. The factors in Eqs. (f6"0"]) and (|6ip . 
except for L and /, represent the mechanical responses of 
a mirror and a bar. The response of a bar is typically 
about 10-times larger. 

The amplitude of the force F n in Eq. (|54[) is propor- 
tional to A. This quantity depends on only the energy 
loss process of exotic particles and the material of the 
mirrors and the bar resonators. This is evaluated as 



.4 



7.3 x 10~ 9 N 



0.42 



Y 



1 - 2(7 



4 x 10 11 Pa/ V5 x 10- 6 /K 
7.9 x 10 2 J/kg/K 
C 



1 dE 



3 GeV/(g cm" 2 ) P dl 



(62) 



In the quantities of Eq. (j6"2"|) , only the linear thermal- 
expansion coefficient, a, strongly depends on the material 
(79] ]. The values of the coefficient a for fused silica and 
sapphire at 300 K, are 5.5 x 10~ 7 /K and 5.0 x 10" 6 JK, 
respectively. The coefficient a of the alloy A15056 [801, 
which is the most popular material of bar resonators @, 
0,&i, is 2.3 x 10~ 5 /K. 

From the above discussion, the advantages of interfer- 
ometers, the higher strain sensitivity and wider observa- 
tion band, are canceled by their larger detector size, be- 
cause exotic-particle detectors must have good displace- 
ment sensitivity, not strain sensitivity. The larger me- 
chanical response (about 10 times) and linear thermal- 
expansion coefficient (several or several tens times) of 
bar resonators enhance the sensitivity. The typical S/N 
of interferometers is obtained from Eqs. ([53"|) , ([58"]) . (|60[) . 
and 



S/N int - 10 1 



5.0 x 10- 6 /K 
1 



1 dE \ f«*\ 
277— 1 ■ (63) 



\3 GeV/(g cm 2 ) P dl 

The S/N of bar resonators is evaluated from Eqs. (|57ll . 
(52), (53, and (S3): 



S/N bax = 3 x 10 2 



2.3 x 10- 5 /K 
1 



ldE\ 
2VT— 1 ■ (64) 



\3 GeV/(g cm^) P dl 

The sensitivity for an exotic particle of bars is a few tens 
or a few hundreds times better than that of interferome- 
ters [8l| . The sensitivity of bars in the above discussion, 
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Eq. (|59|) , is based on the goal sensitivity. The current 
sensitivity is 10-times worse than it [77| . The current bar 
resonators are the better exotic-particle detectors than 
the future second-generation interferometers as LCGT 
HH and Advanced LIGO 

It is difficult to improve the sensitivity of interferom- 
eters for exotic particles. One reason is that in order to 
enhance the signal, the mechanical response and the co- 
efficient of thermal expansion of a mirror must be larger. 
Equation (|60[) implies that a smaller mirror and beam 
yield a larger mechanical response. However, this strat- 
egy enhances the amplitude of the displacement noise, 
and the S/N does not increase. A smaller mirror in- 
creases the radiation-pressure noise (cx R~ 3 ). A smaller 
beam and a larger coefficient of thermal expansion in- 
crease the amplitude of the thermal noise caused by ther- 
moelastic damping in the mirror substrate (cx a/ro 3 ^ 2 ) 
83|. Although mirror cooling reduces the thermal noise 
4a, |46j, l47| , S/N does not become larger, because the ex- 
citation by an exotic particle becomes smaller than that 
at room temperature, as shown in Sec. IIII At 



VII. CONCLUSIONS 

We obtained a general formula for a mirror vibration 
caused by a cosmic-ray particle, and studied the effects 
in typical cases of interferometric experiments. This for- 
mula reveals differences in the responses of resonators 
and interferometers against cosmic-ray particles. In the 
case of resonators, the contribution of the resonant vi- 
bration is dominant. On the contrary, in the case of 
interferometers, the motion of the centers of resonant vi- 
brations must be taken into account. Although the ef- 
fect of cosmic-ray particles of resonators is independent 
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